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SUMMARY 

Two-point , two-time correlation equations are ob- 
tained by considering the Navier-Stolces equations 
f or two points in a turbulent fluid at two different 
imes . By neglecting the triple correlations in the 
equations , a solution is obtained for the space-time 
velocity correlation in the final period of decay 
The analysis is extended to earlier times by con- 
sidering the Navier-Stokes equations at three points 
in the fluid at three different times. The resulting 
set of equations is made determinate by neglecting 
the quadruple correlations in comparison with the 
triple correlations , as in a previous paper by the 
author which considered correlations involving only 
one time. 

The diffusion of particles from a source in a 
decaying turbulent field is calculated approximately 
by assuming that the velocity fluctuations are small. 
The theoretical results are compared with experi- 
ments for diffusion from a line source in a decaying 
turbulent stream. 

INTRODUCTION 

Most of the theoretical work on homogeneous 
turbulence has been based on correlations between 
fluctuating quantities at several points in a fluid 
at a single time (e.g., ref. 1). Correlations in- 
volving several different times as well as several 
points in the fluid are also of considerable interest 
and have been studied by several authors (refs. 
2 to 8). These studies were concerned mostly 
with the kinematics of space-time correlations, 
although some aspects of the dynamical problems 
were also considered. In connection with the 
dynamical problem, Bass (ref. 5) set up the space- 
time equivalents of the K firm h n -H o w ar t h equa- 
tion (ref. 9), but no solutions were obtained. 


This paper is concerned primarily with the 
dynamical problem. First, a solution for the final 
period is obtained b.y neglecting the triple corre- 
lations in the two-point, two-time equation. A 
similar solution was obtained by Batchelor and 
Townsend (ref. 3) by use of a method that con- 
sidered unaveraged velocities rather than the 
two-point, two-time equations considered here. 
However, the method used in this report is more 
convenient for extension to earlier times. The 
extension to earlier decay times, or to higher 
Reynolds numbers, is made by retaining the 
triple correlations. An expression for these corre- 
lations is obtained by neglecting the quadruple 
correlations in a three-point, three-time equation. 
Solutions for still earlier times could be obtained 
by considering the turbulent fluid at a larger 
number of points and times. This procedure is 
analogous to that used previously by the author 
for multipoint correlations at a single time (refs. 
10 and 11). 

By assuming that the turbulent fluctuations are 
sufficiently small for squares and products of the 
fluctuations to be negligible, turbulent diffusion 
from a source is calculated approximately; it can 
be shown that the Lagrangian correlation and the 
Eulerian time correlation are essentially equal for 
this case. The possibility of replacing the La- 
grangian by the time correlation at a point has 
been suggested by Burgers (ref. 4). Recently 
Baldwin (ref. 12) obtained an experimental indi- 
cation that this is a reasonable approximation. 

In the next section the space-time correlations 
for the final period are considered; a higher order 
approximation for earlier times is taken up in later 
sections of the paper. 
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SYMBOLS 

indicates substantial derivative 
energy spectrum function 
arbitrary functions 

constant that depends on initial con- 
ditions 


TWO-POINT, TWO-TIME CORRELATION AND 
SPECTRAL EQUATIONS AND APPLICATION 
TO FINAL PERIOD OF DECAY 

For obtaining the two-point, two-time correla- 
tion equations, first write the Navier '-Stokes 
equations for the points P and P' separated by 
the distance vector r and the time increment 
At: 


integer 

but 

, d (« ,«*) 

_1 d/> b 2 Ui 

(1) 

points 

~dt ~ 


p bii V bx k bx t 

instantaneous pressure 
longitudinal space-time correlation 

bit' 

, yyy) 

1 bp' b 2 u' 

P bx' i +V bx' k b/ k 

(2) 

coefficient, defined by eq. (21) 

bt ^ 

h bx' k 


distance vectors 

dimensionless time, defined by eq. 
(43) 

dimensionless time increment T'—T 
dimensionless time halfway between 
T and T 
times 

time halfway between t and t 9 
reference time 

time increments V — t and t" — f } 
respectively 

instantaneous velocity components 
component of velocity in y-direction 
given by eq. (38) 
space coordinates 

distance in ^-direction that a fluid 
particle originally at y — 0 travels 
during time interval t f ~t x 
space coordinates 

Fourier transforms defined by eqs. 

(27), (28), and (29) 
constant that depends on initial con- 
ditions 

angle between k and k' 
wave number vectors 
equals d^d^d^ 

microscale based on space interval 
microscale based on time interval 
kinematic viscosity 
den si ty 

Fourier transform defined by eq. (7) 
or (S) 


Subscripts: 

L Lagrangian 

iJjkJ tensor subscripts that have values 1, 2, 
or 3 and designate coordinate directions 
Superscripts: 

referring to points P' and P rT 


where the subscripts can take on the values 1, 2, 
3 and a repeated subscript in a term indicates a 
summation. The quantities u* and u'j are instan- 
taneous velocity components, x t is a space 
coordinate, t is the time, p is the density, v is 
the kinematic viscosity, and p is the instantaneous 
pressure. Multiplying the first equation by u'j, 
the second by u if and taking space averages 
result in 


buiu' 

bu,u' i U t _ 

1 bpu’i 

b 2 Ui u' 

bt 

bx k 

p bx t H 

bx t bx k 

buy 

buiii’u't 

1 bp' Uj 

b^y 

bt' 

bx' k 

p w 

H ' bx’tbx'k 


(3) 


(4) 


where the fact that quantities at x * and t are 
independent of j' and t f was used. By introducing 
the transformations b/bxt= — b/br iy b/bxi~b/br 1} 
(bjbt) t ’~ (b/bty&i—b/bAt, and b/bt' — b/bAt } which 
are obtained by writing a correlation as a 
function of r t} t , and At and differentiating, the 
following equations are obtained from (3) and 
(4): 


dUiU; , b 77 7 T 7 , 777 b tt 7 

~dT + dF u J u “ u, (~ r ~ M > t + At )-^jr u t u k u,(r,At,t) 


1 b 


1 c> 


=- v PX,~ ~ pu' i (-r-At,t+At)+2. 

P br t p br } br k br t 

(5) 


bu,u' 


bAt ^-W t {-T,-At,t+M) 


1 b 


b 2 UjU' 


Equations (5) and (G) are the space-time equiva- 
lents of the Karman-Howarth equation. They 
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wore obtained in a slightly different form, for 
the ease of isotropic turbulence, by Bass (ref. 5). 

In order to convert equations (5) and (6) to 
spectral form, the following three-dimensional 
Fourier transforms are introduced: 

u ( Uj(r,M,t)= J Vti (K,M,t)e ik * r d« (7) 

=J *#,(«,**, *)«'* ■ r dK (8) 

pu i {r,M,t.)= J \j(K,At,t)e ,K r d»c (0) 

where #cis a wave number vector and dfc-df^d^d^. 
By introducing these transforms, equations (5) 
and (6) become 


(12) and (13) are neglected, and the following 
solutions are obtained: 

<p tt = ( 14 ) 

(15) 

Tn order for these equations to be consistent, 

E=J(k) C«-'o) ( 1 6) 

where the energy spectrum function E—^vt^tpu 

has been introduced. Evaluate J(k) by letting 

E=JtK A fi* when k is small (Lin, ref. 13). This 
gives 

/(T== J(/ e -^Q-h+l A/) (17) 


^~-{-iK k (pj ki (—K,—M,t+At) — iK k <(i ik j(K,At,t) 

=~iic t \ j — - k,—M , t-\-At) — 2vK 2 <f> ij ( 10 ) 

P P 

§gr*W(-«-*M+A0 

= — I K;X(( — K,— At,t + M) — VKtlpij (11) 


where Jo is a constant that depends on initial 
conditions. For A/ = 0, equation (17) reduces to 
the usual expression for the energy spectrum 
function in the final period, which involves only 
one time. By integrating equation (17) with 
respect to k } the time correlation is obtained as 


UjU'j 

2 


Jo 

32 ( 2 tt ) 1/21 ' 



-'•+5") 


(IS) 


In order to convert the tensor equations (10) and 
(11) to scalar equations, contract the indices i 
and j: 

^f-\-2vK i <pu=iK k lf>,ki(K,M,t) 




and, for isotropic turbulence, the longitudinal 
space-time correlation is 


('-'"+1 a, r’ 


exp\ 


Sr ^ t — ) 


09) 


— Af,t+A<) (13) 

The pressure terms drop out of these scalar equa- 
tions because of the continuity relation 5 ?/*/&/< = 
a?/.;/ajp;=0 and the relation — d/dr'* (see 

eqs. (3) and (4)). 

FINAL PERIOD 

Equations (12) and (13), as they stand, contain 
too many unknowns for solutions to be obtained. 
For the final period of decay, however, the triple 
correlation or inertia terms should be negligible 
compared with the double correlation terms. 
Thus, the terms on the right sides of equations 


Equations (18) and (19) again reduce to the usual 
expressions involving only one time if Af — 0 
(e.g., ref. 1 , p. 94). 

If a new time t m , which lies halfway between t 
and t'(i m =t+At/2) is defined, then At does not 
appear explicitly in equations (17), (18), and (19), 
and U\u[ is a function only of t m and r. For 
instance, the longitudinal space-time correlation 
becomes 

Uitt'i (r,t m ) =4g ( 2°) , /2 v ~ tn 0 ) ~° n 

(20) 
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However, the correlation coefficient is not inde- 
pendent of At. The longitudinal space-time cor- 
relation coefficient R n is defined as 


R w — 


uvu\{r,t m ) 




1/2 


( 21 ) 


and, for the final period, becomes 

(22) 

Batchelor and Townsend (ref. 3) previously ob- 
tained this equation but by a different method, 
which considered unaveraged velocities. A dimen 
sionless plot of 7? n is presented in figure 1. The 
values of 7?n decrease as time interval At increases. 
This is similar to the variation of correlation 
coefficient with distance and would be expected 
physically. The curves go to zero at a finite value 
of At(Af/2^=t m ~to) whereas, as r increases, they 
go to zero only at r=oo, The point where the 
value of Ii n is zero corresponds to the point where 
one of the velocity fluctuations becomes infinite. 


These curves for the final period would not, of 
course, be expected to be accurate in the vicinity 
of that point. The singular behavior could be 
avoided, for positive values of At, by evaluating 
the correlation coefficient at t rather than at t m , 
as in equation (19). However, this coefficient 
would not be symmetric with respect to At. 

A microscale X,, which is based on time interval, 
can be defined by analogy with the usual micro- 
scale X, which is based on space interval. Thus, 
X* might he defined as 

1 = 1 /W? n\ 

X? u*\dAt 2 J 0 

where u 2 =~u i u i for isotropic turbulence. The 
ratio of X 2 to X 2 is then 

/ X t \ 2 il 2 ( d 2 7? i i/dr 2 ) o 

W (^RnldAt 2 ), 

For the final period this becomes 

\\J 1 20(2w) :/2 v 7/2 

Thus, X,/X in the final period is a function of decay 
time as well as of J 0 and v. Calculation of X t /\ 



Figt re 1. Variation of longitudinal space-time double-velocity correlation coefficient in final period (eq. (21)) with 
space and time intervals. Correlation coefficient evaluated at = M/2. 
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from this equation and the experimental data of 
Batchelor and Townsend in the final period (ref. 3) 
indicated values of that ratio on the order of unity. 

THREE-POINT, THREE-TIME CORRELATION AND 
SPECTRAL EQUATIONS 

To obtain the three-point, three-time correlation 
equations, write the Navier-Stokcs equation at the 
points P, P', and P" separated by the distance 
vectors r and r' and the time increments At and 
At'. The vector configuration is shown in figure 2. 



three-time correlation equations. 


The first two equations arc the same as equations 
(1) and (2) with the subscripts k replaced by l. 
The third equation is 


du'*' 

dt" 


d , „ 1 bp” 

dr," (m * ‘ ) pbx k ' 


dr," dr," 


(23) 


By mul tipping the first equation by the 

second by UiU k , the third by utfij and by taking 
space averages, the following correlation equations 
can be constructed : 


77 U i U i U k - UiUjU* Wi — 377 UtUjUt Ui 

Ot OVi Ot'i 


+ ~ uMuWl+^ r 7 Uiulul'ul' 


1 / d / // d / // 

=-A-W i 71UiWk -M pUiU * 


+ 2p 


( 


b 2 u,u'u k bhijU^l 


dr, do 


dr, do 


yufUjUh \ , 0 „ 

dr'M ) 


^ + UiU'UtU’, 

1 a rr 


1 v / f7 i a WiUjUi- 

^—7> «««*+«• ^ >. (25) 

p dr/ do dr, 


^7 UtUfUt u, 


1 bp" Uiii’i , d 2 u,u,V*' /0 ^ 

= YV - > - (20) 

p C)/’* O/'jOf ^ 

where the following transformations were used: 


d d 5 d c) _d^ 

dx t dr t dr/ dx'^ dr/ dxY dr/ 


(d<X,,.. (df)a..A,' dAf bAt'’ 

b b b _ d 
bt’~~bAt’ bt " bAt' 


Equations (24), (25), and (26) are the three-point, 
three-time correlation equations. In order to 
convert these equations to spectral form, the 
following six-dimensional Fourier transforms can 
be defined: 


UiU'jUtir.Atyr'At'.t) 

O ” p ijk (K,At ,k' ,A t' ,t)e HK • r+ *' ■ r0 dKdK' 

- 03 

(27) 


UiUtUjUt (r,At,r' ,At' ,t) 

= f’ r ^ad^AtyAt',^"^'"'^^' 

C/ — CO 1/ — CO 

(28) 


Oa 


a Jk e t(K ' T+ *' • r ' , d«dK' (29) 


Also, 




= Ujupi'iU k ( — r, — A t,r' — r,At' — A/,£+Af) 

= f ^jiiic(—K — K', — M,K',At'—At, 

%J — CO 

‘ r+ * * r (30) 


and 


r',— AP,r— r' , At— Ap,t+AP) 

= f K— If',— At',K,At— At', 

<+A#')e 1( “ r+lt '' r ' ) dKdK' (31) 
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Similar expressions can be obtained for pressure 
correlations. 

Substituting tlie preceding relations into equa- 
tions (24), (25), and (26) gives, for the three- 
point, three-time spectral equations: 

A P tik +2v(K i +K,Kl + K n )P, j t=i(K l + K l )P i , Jk 
-iK l Pnd-K-K'-At,K , > At'-At ) t+ At) 
—iK'i@kui(—K—K',—At',K,At—At',t-\-At') 

~p f~ + K — K' , — At ,K , 

At' — A/,£-|-A0 + ?K* a i7( — k—k' } — At' } k, 

At-At',t+At')] (32) 

i'Klftjiiki. K K , At } K , 
At'-At,t+At)-- ?K,a,*( — k—k', 

P 

—At.K'At'—Mt+At) (33) 
^f7+ VK ' 2 PtJX = — X— *’>— &t', K, 


At-At',t+At')—-iKla u (—K-K',—At',K, 

P 

At-At',t+At') (34) 


SOLUTION FOR TIMES BEFORE FINAL PERIOD 

For (lie final period of decay, a solution was 
obtained by neglecting the triple correlations in 
the two-point, two-time equations. Similarly, a 
solution applicable to times before the final period 
could be obtained by retaining the triple correla- 
tions and neglecting the quadruple correlations in 
the three-point, three- time equations. A fuller 
discussion of this procedure is given in references 
10 and 11. In reference 10 it is shown that, if 
terms corresponding to quadruple correlations are 
neglected, the terms corresponding to pressure 
correlations must also be neglected. Thus, if all 
the terms on the right sides of equations (32), (33), 
and (34) are neglected, the equations can be 
integrated between t 0 and t to give 

Pijk=juk(*<y At At') exp [— 2F(« 2 +fc,K i '-|-#c /2 )(^— 1 0 )] 

Pok=gijk(K,*’ ,t A*') exp (~vk 2 AI) 

PijK=hijk(Ky ,t ,At) exp (— vk' 2 M ') 


In order for these equations to be consistent, 
KkPiik— KklPiik(K,K')]o expj^ —2v (t — t 0 +~Atj 

+ KK'(t-l o) cos e+K ' 2 (i-<„+^Af/)J (35) 

when 1 the indices i and j have been contracted 
and the equation has been inner-multiplied by K k 
in order to convert it to scalar form. The sub- 
script 0 refers to the values of 0 iik at t~t &i At— At' 
= 0; and 6 is the angle bet ween k and k . 

In order to connect (3 tik and <p iki) let r' = At' = 0 
in equation (27) and compare with (8) to obtain 


I. 


<Pi hi (K,At,t)= fi iik (K t At,K ,0,0d«' (36) 


Substitution of equations (35) and (36) into (12) 
results in 

^+2vSE=W (37) 

where 

IZ=2wk 2 <pu 

and 


r* co 

TC=J / (/3 IiJt ) o ( 2 tt) 2 A' 2 exp{—2v\n i (t — t l) 
+At/2) + K n (t — t il )]} ■ er]> [— 2vkk'(1— Q 


s 0](l(eos0)J<lK'-|~J^ [i( «> 

«')] o(2x) W 2 erp-f —2v (<— <o+y) 


+K ,2 (t-t 0 +At ) 


]}■[/: 


exp [—2vKK'(t — t 0 


+A0]d(cos 6 ) J d k' (38) 

where d k' is written as — 27r/ 2 d(cos 0)d*'. The 
quantities (/S iU ) 0 depend on the initial conditions 
of the turbulence. In order that these results 
will reduce to previous results for Af=0 (ref. 10), 
let 

(2?r) 2 ?**&<*(«, k')o=-- |/3 0 (kV 6 — A' 4 ) (39) 

Then, 

(2ir) 2 -«')u=-^o(A /6 -«V 4 ) (40) 
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Substituting equations (39) and (40) into (38) and 
carrying out the integrations with respect to 0 and 
k result in 


W= 




c 


105k 15 . 45 k" 

/ ,i {t-t 0 Y Ji+ p Vit (t-t 0 y /2 


19k 1 


3k 12 "| 

> / 2 («-*o ) s/2 S'V-t o) 3/2 J 


Sim & exp 

256 2 P 


t—t 0+- - At 


[ 


105k 6 


y vli (t-tu+Aty /i+ v f,2 (t-t 0 +Aty n 
19k 10 3k 12 


45k 8 


«] (41) 


v m {t-t Q +MY r% v V2 (t-t 0 +Aty' 2 _ 

Equation (41) reduces to the expression for the 
energy-transfer function involving only one time 

(ref. 10) if At= 0. Note that J Trd«=0 only 

for At=0 . 

Substituting equation (41) in (37) and inte- 
grating with respect to t result in 


E= 


37T 


[ 


exp^- 2vJ (<-<o+|a<)J 

[-§-’( 

12 v 2k 8 




ft 


exp 


256r 2 
15 v2k 6 

v V2 ( t -to) m v 5f2 (t-to ) 5/2 


t — t Q -{ -- At 


)] 


IGv^k 12 


'3v 3/5 (<-< 0 ) 3/2 ’ + "3»’ 1/li («-«o) ,/ * 

32k 13 ,,/ r(<-t 0 ) ,/2 '' 


01 


„-I/2 


ft 

~256r 2 


[■ 


erp — | k 2 v (t— < q+" A <) 


15v2k 6 12 v 2k 8 

V /2 (t-< 0 +Af) 7/2 r 5/2 (f-t 0 +Af) 5/2 


where 

F(co)=e _ “ J J* e J \\x 

„{«]’ or .[ »a-<.+A0 ]- 

For evaluating the function of integration (with 
respect to time) in equation (42), the theory of 
Lin (ref. 13) or of Batchelor (ref. 1) is used. 
According to those theories the coefficient of the 
first term in the expansion of E in powers of k is 
independent of time, whereas the other terms may 
not he. Thus, the function of integration, since 
it cannot be a function of time, is set equal to the 
first term of the expansion, which is J 0 k*/3tt. The 
theories of Tan and Batchelor are based on the 
assumption that correlations arc exponentially 
small for large values of r. This is consistent 
with the results of the present theory (refs. 10 
and 11) and does not seem to be inconsistent with 
some results of Batchelor and Proudman (ref. 14), 
if the effects of the singularities arising in their 
analysis are assumed to be negligible. According 
to Batchelor and Proudman, correlations, in 
general, would be expected to be negative power 
functions of r for large values of r, and all terms in 
the expansion of J? would, in general, be functions 
of time. However, their results by no means 
rule out the possibility of exponentially small 
correlations or of a constant first term in the ex- 
pansion of E. The present theory gives ex- 
ponentially small correlations for large r (refs. 10 
and 11) and is consistent with a constant first 
term in the expansion of E for evaluating the 
function of integration in equation (42). 

By using the relation 


UiUi(At } t) 

2 



E<\k 


7^ k 10 16 v ^k 12 

+3 *» /2 (< — t 0 +Aty'* + 3r ,/2 « - t 0 +At) m 


32k 13 H'tt-to+AQ 

t'vL 2 



(42) 


and introducing the dimensionless time 

/r 


(43) 


5S95G1 — 61 2 
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(ho dimoRsioiilcss time correlation becomes 
ffr 1 (t-l.- ,\r\ in — 

jw9„i/9 2 32(27r) l/2 \ J 512 

25/(3 y/3) 25/(3y3) 


T 7/2 (7 , +qA7 , y /2 (T+AT 1 ) 7 ' 2 (t+’jAT 1 ) 
140/(9^3) 140/(9 V '3) 


/ 9 \ 9/2 

7W2/T+|a7M 

(7 , + AF) 5/2 (r+lA7’) 9 ' 2 

, 245/(27 y3) 

, 245/(27 \3) 

~ / 9 \lI/2 

T 3,i (T+^AT) 

r /I \ n/2 

(T +A7 T ) W2 (7 , + ^ AT 7 ) 

. 61 60/ (sly's) 

1 6160/(81y/3) 


( O \ 13/2 1 / 1 \ 13/2 

AT) (T+&Ty f *(T+± AT) 

32 » 1 . 3 * 5 . . . ll+2n 

O 


3 (2n— 1) (ft— T) ! 2 3n+13 


2«- 1 
i 2 


(T+A? 7 ) 


2n+ 13 


(t+IaT’) 2 ( t +^ at ) 


h (44) 


where 0! = 1. For A7 T =0 J this becomes 


ffi 9 UiU t (t) 

2 


7 1 -5/2 

55(SF+ 0 - 2m7 ’" 


( 45 ) 


which is the same as the expression for the turbu- 
lent energy obtained previously (ref. 10). As in 
the case of the final period, an average time t m = t 
+ A?/2 can be introduced. Equation (44) becomes, 
when written in terms of the dimensionless T m 
rather than T } 

UiUiiAfJ 1 r _ 5/2 

Ji 4/s V /9 2 32(2 ir) m m 


7T 

AT2 


25/(3 V3) 


25/(3 V3) 

(r.+lir)' ,, (r.-lA2-) 


7/2 


• Y (40) 


This expression for UiU\ does not become independ- 
ent of AT 1 when written in terms of T m) as was 
the case for the final period. However, the 



Figure 3. — Variation of time double-velocity correlation coefficient (eq. (47)) with dimensionless time interval and 
decay time, and comparison with values for final period. T defined by equation (43) ; correlation coefficient evaluated 
at T m ~ T-\- AT/2, 
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expression is still symmetric with respect to AT. 
This is in agreement with the results of Mcecham 
(ref. 7), who came to that conclusion by kinemat- 
ical considerations. 

A more physically meaningful quantity than 
mUi is the time correlation coefficient, defined as 
follows : 


Ji it = r 




h U f (<m+2 U ‘ U I ( < ™ — | A *) 


(47) 


A plot of Iiu, obtained from equations (45), 
(46), and (47), is shown in figure 3. The final- 
period contributions are shown as dashed curves 


for comparison. As in the case of the final period, 
the values of Tin decrease with dimensionless time 
separation. That is, of course, the type of be- 
ll avior that would be expected on physical grounds ; 
the time correlation coefficient can be considered 
as a measure of the sameness of the velocities at 
different times at a point in much the same way 
that the space correlation coefficient provides a 
similar measure for velocities at different poin ts at 
one time. It is of interest that the values of Uiu'i 
by themselves do not exhibit this behavior; in 
fact, they increase rather than decrease with time 
separation, as shown in figure 4. This unusual 
variation is apparently due to the nonlinear decay 
of the turbulence with time and would not be 
observed for stationary turbulence. For the 



Figure 4.— Variations of dimensionless double-velocity correlation with time interval and decay time. Correlation 

evaluated at T m . 
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decaying case, it appears that the correlation 
coefficient, as defined b y eq uation (47), is much 
more meaningful than u f uj by itself. The dis- 
cussion of the time microscale given in connection 
with the final period would, of course, hav e been 
meaningless if it had been based on UtUt rather 
than on the correlation coefficient. 

Comparison of the dashed with the solid curves 
in figure 3 indicates that the general effect of the 
higher order inertia terms in the correlation 
equations is to decrease the correlation coefficient 
at a given value of time separation. This is 
opposite to the corresponding effect for space 
correlation (see, for instance, the experimental 
results of Stewart and Townsend (ref. 15) or the 
theoretical results of Deissler (ref. 11)). It is 
possible that the reduction of correlation coefficient 
by inertia terms is caused by the nonhomogeneity 
of the turbulence with time. 

APPROXIMATE CALCULATION OF TURBULENT 
DIFFUSION FROM A SOURCE FOR SMALL 
VELOCITY FLUCTUATIONS 


Similarly, the Lagrangian correlation is expanded 
as 


[K<MOk=*' 2 (<) + KO t '_ t (C-0 

(49) 


The substantial or particle derivative can be 
written as 


Dv(n_Mn , „xn wn , ,xo 

Dt' dt' + d* +l by + dz 
For small velocity fluctuations, 


(50) 


Dv(t') ^dv(t') 

Dt ' ® dt' 

Also 

n*v(t') b Dv(t') c Mt') 
~ Dt' 2 =bt' Dt' = dt' 2 


(51) 


(52) 


The time correlations considered in the pre- 
ceding sections were concerned with velocities at 
different times at a fixed point in the fluid 
(Eulerian correlations). On the other hand, 
calculating the turbulent diffusion of particles 
from a source usually involves the Lagrangian 
correlations, which are based on the velocity of a 
moving fluid particle at different times, rather 
than on the velocity at a point. For small 
velocity fluctuations, however, it has been sug- 
gested by Burgers (ref. 4) that the two correlations 
should not differ greatly. This can be shown as 
follows : 

First, consider the Eulerian time correlation 
v(t)v(t') } where v is the component of the velocity 
in the ^/-direction ; similar results could be obtained 
for the other velocity components. The Eulerian 
correlation can be expanded in a series as 




=v>(t)+m 


+ \ v{t) \jw^\ t .J t '~ ty+ • ■ • (48) 


and so on for higher order derivatives. From 
equations (48), (49), (51), and (52) 


v(t)v(t')^[v(t)v(t')] L (53) 

is obtained, which was the relation to be proved. 
It should be noted that relation (53) is most 
accurate for small values of t' — t as well as for 
small velocity fluctuations, inasmuch as t he approx- 
imate relation (51) had to be applied a greater 
number of times to the higher order derivatives 
in equation (49) than to the lower order ones (sec 
eq. (52)). 

It should also be emphasized that equation (53) 
was obtained for the case of no mean motion. 
Thus, Eulerian time correlations measured with 
a stationary instrument in a moving stream will 
probably differ considerably from the Lagrangian 
correlations. However, if the instrument is mov- 
ing with the stream, the two correlations will be 
approximately equal if the turbulence level is not 
too high (see ref. 12). 

Next, the equation for the turbulent diffusion 
of particles originally concentrated at a source is 
considered. The theory of turbulent diffusion was 
originated by G. I. Taylor in 1922 (ref. 16) and has 
since been studied by a number of authors (c.g., 
refs. 17 and 18). The distance in the ^-direction 
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that a fluid particle originally at ?/= 0 travels 
during the time interval t'—ti is 

Y{t')= f *' v(t)dt (54) 

A 

Multiplication of this equation by v(t') gives 


where T x is again the dimensionless time at which 
diffusion begins and AT 2 is the dimensionless time 
during which diffusion takes place (T 2 —Ti). For 
large diffusion times, 


^V 7/9 x^_ 1 

2 d\2irTY 2 


(59) 


v(t')Y(t')=l v(t)v(t')dt (55) 

Taking the particle average over all the “marked” 
particles that were originally concentrated at a 
source at y— 0 and integrating with respect to t' 
result in 

71=2 f 2 f " [v(t)v(t')] L dtdt' (50) 

Equation (50) gives the mean square of the dis- 
tance that the marked fluid particles concentrated 
at y ~ 0 at time fi have traveled by time t 2 . It 
is evidently applicable to decaying as well as to 
stationary turbulence. Note that the double 
integral in equation (50) cannot be converted to 
a single integral as in the case of stationary 
turbulence. 

If the approximate relation (53) is introduced 
and it is noted th at, for isotropic turbulence, 
v(t)i'(t r ) = [u i (t)u i (t r )]/3 ) then equation (56) can be 
written in dimensionless form for isotropic turbu- 
lence as 


^v 7/9 ^_2 r T > r pruMuw) , irrArr , 

^3^ J t x A'V 79 


(57) 


where the dimensionless time T is defined by 
equation (43) and the time correlation is obtained 
from equation (18) or (44) by remembering that 
At=t'—t or AT=T'—T. 


FINAL PERIOD 


For diffusion in the final period of decay, equa- 
tion (57) can be integrated to give 


^V 779 J_ 
JJ° 79 2 9V? 


VZ [(A+at^+tt 2 ] 


vz 


( T '+i AT >y 


That is, the turbulent diffusion distance reaches a 
constant value and becomes independent of AT 2 
for large diffusion times. This differs from the 
case of stationary turbulence, where Y\ increases 
linearly with AT 2 for large diffusion times. The 
reason it reaches a constant value for decaying 
turbulence is that for large times the turbulence 
goes to zero, so that no more turbulent diffusion 
can take place. 

Figure 5 shows dimensionless root -mean-square 
diffusion distance for the final period plotted 
against diffusion time for various values of T l} the 
time at which diffusion begins. The curves have 
considerable curvature at early times but approach 
a linear form for large values of 7V For early 
times, the diffusion distances are much larger than 
those for later times because of the higher turbu- 
lence level at early times. 

TIMES BEFORE FINAL PERIOD 

It might be argued that, strictly speaking, the 
approximate relation (53) should be used only in 
the final period, inasmuch as inertia terms were 
neglected in obtaining it. Inasmuch as no experi- 
mental diffusion data exist for the final period, 
however, some sort of approximation must be 
made for earlier times in order to compare the 
theory with experiment. The results might still 
be applicable for small times of diffusion; experi- 
mental data of Baldwin (ref. 12) for diffusion in a 
fully developed pipe flow indicated that equation 
(53 applies reasonably well for that case, although 
the turbulence probably did not correspond to that 
in a final period of decay. 

For times earlier than those corresponding to 
the final period, equation (44) is used in (57) with 
AT replaced by T f — T. In this case the integra- 
tion was carried out numerically on high-speed 
computing machinery. The resulting plot is 
shown in figure G, where the final-period contribu- 
tions are shown as dashed curves for comparison. 
The higher order inertia terms have a noticeable 
effect on the diffusion at early times; at later 
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Figure 5. — Predicted root-mean -square turbulent diffusion distances for final period as a function of dimension- 
less diffusion time and decay time. 


times the effect of those terms becomes negligible, 
and the solid curves approach those for the final 
period. 

These curves should apply to the calculation of 
the width of the diffusion wake from a line source 
in a moving stream. In this case At would be 
replaced by the distance downstream from the 
source divided by the velocity of the mean stream. 
Comparison of the curves in figures 5 and 6 with 
those obtained experimentally (e.g., ref. 19) does, 
in fact, indicate a marked similarity. In order to 
obtain a more quantitative comparison, the 
constants J 0) ft, and t 0 , which depend on initial 
conditions, were evaluated from the decay data of 
Ubcroi and Corrsin and equation (45). Equation 
(45) was found to represent the decay data closely 
when J Q = 1.05X10- 9 ft 7 /sec 2 , ft=1.81X10- 33 
ft I8 /sec 3 , and ^=—0.407 sec. With these values 
for the constants, diffusion data for an early and 
a late time are plotted in figure 7. Included for 
comparison are analytical results for the same 


values of T v . The agreement between theory and 
experiment seems to be good for large values of T x 
and small values of AT 2 , whereas some deviation 
is indicated for other conditions. This might have 
been expected from the nature of the approxima- 
tions made in obtaining equation (53), which was 
used in the analysis. As discussed previously, 
that relation is most accurate for small velocity 
fluctuations (large T } ) and for small diffusion 
times. 

CONCLUSIONS 

The time correlation coefficient in a decaying 
homogeneous turbulent field, when evaluated at 
a time halfway between the limes at which the 
two velocities are considered, decreased with time 
interval in much the same way that space correla- 
tion coefficients decrease with space interval. The 
time correlat ions by themselves, on the other hand, 
were independent of time separation in the final 
period and increased with time separation at 
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Figure 6. — Predicted root-inean-square diffusion distances for times before final period and comparison with final -period 

contributions. 


earlier times, although they were symmetric 
with respect to At. The correlation coefficient 
(eq. (47)) appears to be a much more physically 
meaningful quantity than the correlation for a 
decaying turbulent field. The effect of the higher 
order inertia terms in the correlation equations 
for times before the final period was to reduce the 
value of the correlation coefficient at a given time 
interval below that for the final period. The 
ratio of time microscale to space microscale in the 
final period was a function of decay time and of 
initial conditions. 

By assuming that the velocity fluctuations are 
sufficiently small for squares and products of 
velocities to be negligible, it can be shown that 
the Eulerian time correlation is approximately 


equal to the Lagrangian correlation. Turbulent 
root-mean-square diffusion distances were calcu- 
lated by using this approximation and the equa- 
tions for the time correlation obtained herein. 
The agreement between theory and experiment 
was good for large decay times (low turbulence 
levels) and for small diffusion times; for other 
conditions, some deviation was indicated. This 
was apparently due to the assumption of the 
equality of Eulerian and Lagrangian correlations, 
that assumption being most accurate for small 
velocity fluctuations and short diffusion times. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 16 , 1060 
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Figtre 7. — Comparison of theory and experiment for width of turbulent diffusion wake from a line source in a decaying 

turbulent stream. 
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